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Image  Charge  Forces  Inside  Conducting 

Boundaries 

Mark  D.  Tinkle*  and  S.  E.  Barlow* 


* Dynamics  Technology,  Inc.,  21311  Hawthorne  Blvd.,  Suite  300  Torrance,  CA  90503-5610 
t W.R.  Wiley  Environmental  Molecular  Sciences  Lab,  Pacific  Northwest  National 
Laboratory  P.O.  Box  999  (K8-88)  Richland,  WA  99352 

Abstract.  Evaluation  of  the  force  on  a charge  particle  due  to  the  surface  charge  distri- 
butions it  induces  in  nearby  conductors  is  generally  a challenging  problem.  We  show  these 
forces  can  be  evaluated  explicitly  and  fairly  simply  in  a number  of  elementary,  but  important 
cases. 


INTRODUCTION 

The  force  on  a charged  particle  due  to  the  charges  it  induces  on  nearby  conducting 
surfaces  is  a factor  in  many  scientific  measurements.  The  effects  of  these  forces  are 
particularly  important  for  non-neutral  plasmas.  Until  now,  general  techniques  for 
quantifying  image  charge  forces  have  not  been  available.  An  important  advance 
was  made  in  a recent  paper  by  Fine  and  Driscoll  [1]  which  addressed  the  lowest 
order  solutions  to  this  problem  for  infinite  cylindrical  geometry  including  finite 
charge  length  effects.  Here,  we  generalize  those  results  considerably. 


IMAGE  CHARGE  PSEUDOPOTENTIAL 

The  measurable  physical  variable  associated  with  induced  surface  charge  is  the 
force  it  produces  on  the  point  charge.  We  would  like  to  express  this  force  as  the 
gradient  of  a function  of  the  point  charge’s  position.  There  is  a subtlety  involved 
in  this  potential  function  that  is  best  illustrated  by  a simple  example  using  the 
method  of  images. 


Motivating  example 

If  a point  charge  of  strength  q is  at  a position  ( x,y,z ) = (a;o,0,0),  (x0  > 0)  with  a 
grounded  conductor  in  the  y-z  plane,  we  know  from  the  method  of  images  [2]  that 
the  electrostatic  potential  for  x > 0 can  be  expressed  as  the  sum  of  the  potentials 
from  the  point  charge  and  a fictitious  image  charge  of  strength  —q  located  at 
(x,y,z)  - (—a:o,0,0): 
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47rC°  \\/[(z-*o)2+2/2  + 22)]  y[(x  + .T0)2  + y2  + Z2)]J 

for  x > 0.  This  is  the  potential  that  determines  the  force  on  an  infinitesimal  point 
charge,  which  is  to  say  a point  charge  whose  own  induced  surface  charges  produce 
negligible  forces.  Can  we  use  it  to  find  the  force  on  the  charge  q ? We  must  first 
drop  the  potential  produced  by  the  charge,  since  it  cannot  exert  a force  on  itself. 
This  leaves  the  potential  of  the  image  charge,  which  evaluated  along  the  x axis  is 
just 

*/(*, 0,0)  = (2) 

47T€g  X "f*  Xq 

for  x > 0.  Two  courses  of  action  present  themselves  for  finding  the  x-component  of 
the  force:  take  d$i/dx  and  evaluate  at  x = x0,  or  evaluate  at  i = x0  and  take 
d/dx0  of  the  result.  The  results  differ  by  a factor  of  two: 
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dx 

X=X  o 
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dx0 

47tc0  2xq 

The  force  can  be  found  directly  from  Coulomb’s  Law  using  the  image  charge[2]: 

,5) 

so  at  least  in  this  case,  we  have  a choice  of  formulae  for  the  force  in  terms  of  the 
potential: 

F(x0)=-q^-  x (6) 

£=ro 

or 

pi  , _ -q9[^/Ur,] . m 

F(  o)  2 dx0  () 

The  first  choice  seems  like  the  natural  one,  but  it  requires  us  to  keep  a function 
of  two  position  coordinates  (x  and  x0).  The  second  choice  contains  an  unfamiliar 
factor  of  1/2,  but  the  force  is  reconstructed  from  a simple  function  of  one  position 
coordinate  (x0).  This  is  the  useful  form,  which  we  will  now  derive  more  generally. 


General  derivation 

Next  we  find  a general  expression  for  the  force  on  a point  charge  q in  the  presence 
of  both  fixed  charges  and  conductors  at  fixed  potentials. 
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where  <f>  is  the  ordinary  electrostatic  potential  due  to  the  fixed  charges  and  poten- 
tials, and  $ is  the  potential  due  to  charges  induced  by  the  point  charge,  evaluated 
at  the  point  charge  coordinate. 

From  the  definition  of  the  electrostatic  potential  <j>,  the  work  required  to  bring  a 
point  charge  q in  to  position  x from  infinity  (where  <f>  = 0)  while  holding  the  charges 
producing  <j>  fixed  in  place  is  W — q<j){x).  The  force  on  the  particle  is  the  negative 
of  the  gradient  with  respect  to  its  position  coordinate  of  the  total  system  energy 
under  the  circumstances  considered.  If  the  charges  producing  4>{x)  are  indeed  fixed 
and  there  are  no  other  energy  terms  involving  the  position  of  the  point  charge, 
then  the  force  is  F — —qS/ej)-  This  is  the  familiar  result  for  an  infinitesimal  test 
charge,  but  it  is  a special  case.  It  is  not  the  definition  of  (f>. 

When  a charge  moves  toward  a conducting  surface,  the  induced  surface  charges 
must  bunch  together,  against  their  mutual  repulsion,  to  maintain  an  equipotential. 
This  adds  a term  to  the  total  electrostatic  energy  of  the  system  that  acts  to  reduce 
the  attractive  force  felt  by  the  point  charge,  by  exactly  a factor  of  two.  This  is 
easily  derived  from  the  well-known  expression  for  the  total  electrostatic  energy 
(excluding  self-energy)  of  n discrete  charges: 

W = l±q^,  (9) 

z j=i 

where  3>j  is  the  potential  at  charge  q3  due  to  all  the  other  charges. [3]  If  j = 1 denotes 
the  point  charge,  and  the  others  are  the  charges  induced  on  a grounded  conductor, 
then  is  zero  for  all  j ^ 1 (the  induced  charges)  regardless  of  the  position  of  the 
point  charge,  and  the  total  electrostatic  energy  is 

W=\q *,  (10) 

where  $ = 4>i  is  the  potential  at  the  point  charge  due  to  theinduced  surface  charges. 
The  force  on  the  point  charge  is  thus 

F = (II) 

where  the  gradient  is  taken  with  respect  to  the  position  of  the  point  charge.  If  the 
conductor  is  at  some  potential  other  than  ground,  then  we  must  add  the  familiar 
term  —qV<j>  to  this  thereby  recovering  eq.(8). 

We  refer  to  $(x)  as  the  image  charge  pseudopotential,  to  emphasize  its  dif- 
ferences from  ordinary  potentials.  Unlike  <f>{x),  T>(x)  is  not  simply  a solution  of 
Laplace’s  or  Poisson’s  equation,  because  its  source  term  (the  induced  charge  distri- 
bution) is  a complicated  functional  of  x,  the  point  charge  location.  Further,  $(x) 
is  only  meaningful  at  the  particular  location  of  the  charge.  As  we  describe  below, 
$(a?)  is  proportional  to  the  nonsingular  part  of  G(x,x'),  evaluated  at,  x = x'.  Thus, 
$(£)  may  be  extracted  from  a family  of  solutions  to  Laplace’s  equation,  but  is  not 
itself  a solution.  Numerous  methods,  include  the  classic  “Method  of  Images”  for 
solving  these  problems  are  discussed  in  Tinkle  and  Barlow. [4] 
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RESULTS 


Generic  limiting  form 


As  a smooth  region  of  a conductor  is  approached  closely,  the  pseudopotential  ap- 
proaches that  of  an  infinite  flat  plate.  From  the  method  of  images.  $ ->  —q/AntuId 
as  d — > 0,  where  d is  the  distance  to  the  surface.  Thus  if  60  is  the  distance  to  a 
conductor  along  the  6 axis  of  a rectilinear  coordinate  system  and  the  surface  is 
perpendicular  to  the  6 axis  where  they  intersect, 


$ -*  - 


<?  1 

4^06,2(1—0’ 


(12) 


as  6 = 6/6)  — ^ 1.  If  the  geometry  is  symmetric  in  6,  we  can  account  for  the  limiting 
behaviors  at  both  walls  with  the  form 


$ ->  - 


g 1 

47rc06o  1 - 0 ’ 


(13) 


as  |6|  — y 1.  The  complete  solution  for  $ along  the  6 axis  will  be  the  sum  of  this 
divergent  term  and  a finite  term.  Expressed  as  a power  series  expansion  about  the 
origin,  it  will  have  the  form 


4>  = 


4™o6o  iS 


(14) 


with  C„  = 0 for  odd  n,  due  to  symmetry.  The  divergent  term  by  itself  lias  C„  = 1 
for  all  even  n.  The  results  we  obtain  for  specific  geometries  can  all  be  expressed  in 
this  form,  with  different  values  for  the  coefficients.  Tinkle  and  Barlow[4]  present 
solutions  of  this  form  for  a number  of  important  cases  including  the  spheical  shell, 
parallel  plates  and  the  rectangular  box. 


Cylindrical  Shell 


By  way  of  example,  we  show  here  the  results  for  a cylindrical  box.  The  calcula- 
tions are  rather  complex, [4]  however,  solutions  of  the  form 


2m  -2 n 


(is) 


can  be  found.  Here  a is  the  ratio  of  the  cylinder’s  length  to  diameter-the  “aspect 
ratio. ’’The  limiting  values  listed  in  Table  1 are  in  line  with  our  generic  discussion 
above.  These  coefficient  terms  are  shown  plotted  in  Figure  1 . 
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Aspect  Ratio  a(Zo/  ro) 

FIGURE  1.  Plots  of  the  the  first  six  terms  of  the  power  series  expansion  of  the  image  charge 
pseudopotential  for  a cylindrical  box  as  functions  of  the  aspect  ratio.  (A)  The  pure  “axial”  terms- 
powers  of  z only.  (B)  The  pure  “radial”  terms-powers  of  r only.  (C)  The  crossed  “radial”  and 
“axial”  terms-powers  of  r and  z. 
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TABLE  1.  Asymptotic  limits  of  image  charge  coefficients  for 
cylindrical  box,  see  text  for  explanation.  *agrees  with  results  of 
Fine  and  Driscoll  [ref  1,  eq(32)j 


a — y 0 
parallel 
plate 

Ath  Order 
Exp’l.  Trap 

a=  1 

a -4  0 
long 
tube 

Co,2 

1.0518 

0.7755 

0.5880 

■ ■ B 

Co, 4 

1.0045 

1.1669 

1.1845 

0 

1.0005 

1.0180 

1.1246 

C2,0 

0 

0.7534 

0.8792 

1.0027* 

cy  o 

0 

11243 

1.0884 

1.0009 

Cs,a 

0 

1 .0572 

1.0020 

1.0003 

Cl, -2 

0 

-0.9157 

-1.0309 

0 

C 4,2 

0 

-0.3470 

-0.4359 

0 

<?2,4 

0 

0.1822 

0.1884 

0 
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